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ABSTRACT 

We use a spherical non-LTE fully line blanketed model atmosphere to fit the 
full multi-wavelength spectrum, including the extreme ultraviolet (EUV) continuum 
observed by the Extreme Ultraviolet Explorer, of the B2 II star e CMa. The available 
spectrophotometry of e CMa from 350 A to 25 /^m is best fit with model parameters 
Tes = 21750 K, log (7 = 3.5, and an angular diameter of 0.77 mas. Our best fit model 
predicts a hydrogen ionizing flux, qq, of 1.59 x 10^^ photons cm~^ s~^ at the star's 
surface and 2290 photons cm^^ s^^ at the surface of the Local Cloud. 

The close agreement between the model and the measured EUV flux from e CMa 
is a result of the higher temperatures at the formation depths of the H I and He I 
Lyman continua compared to other models. The realistic model treatment of early 
B giants with spherical geometry and NLTE metal line blanketing results in the 
prediction of significantly larger EUV fiuxes compared with plane-parallel models. 
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We find that our metal line blanketed spherical models show significantly warmer 
temperature structures, 1-3 kK at the formation depth of the Lyman continua, and 
predict stronger EUV fluxes, up to a factor of 5 in the H I Lyman continuum, compared 
with plane-parallel atmospheres that have identical model parameters. In contrast, we 
find spherical and plane-parallel models that do not include metal line blanketing are 
nearly identical. Our T^s = 21000 K, \ogg = 3.2, spherical NLTE model predicts more 
than twice as many hydrogen ionizing photons and over 200 times more neutral helium 
ionizing photons than a standard hydrostatic plane-parallel LTE model with the same 
stellar parameters. 

Our synthetic spectra are in reasonably good agreement with observed continuum 
and line fluxes from echelle spectra obtained with the Goddard High Resolution 
Spectrograph. While we find agreement between the absolute UV flux of e CMa as 
measured by GHRS and our model atmosphere, these fluxes are ~30% higher in the 
UV than those measured by lUE, OAO-2, and TD-1, in excess of the published errors 
in the absolute calibration of these data. 



Subject headings: stars: individual (e CMa, /? CMa, a Vir), early- type, atmospheres, 
fundamental parameters (radii, temperature, gravity) 



1. Introduction 



Spectroscopy of the bright B2 II star e CMa (HD 52089, HR 2618) dCassineUi et al. 1995| ) 
below the Lyman edge by the Extreme Ultraviolet Explorer (EUVE) satellite is possible due to the 
extremely low neutral hydrogen column density toward this star. As a result of its location at a 
HIPPARCOS distance of 132^g^ pc along an exceptionally rarefied interstellar tunnel extending 
out from the Local Bubble ( Welsh 198S| ), e CMa is attenuated by a neutral hydrogen column 



density of less than 5 x 10^^ cm ^ ( |Gry et al. 1995 ). Consequently, e CMa is an extremely 



important star for its contribution to the local interstellar hydrogen ionization, producing more 
hydrogen ionizing flux than all nearby stars combined ( [Vallerga &: Welsh 1996 ). The large EUV 



flux from e CMa affects the ionization state of the Local Cloud, the region of neutral hydrogen 



concentrated within a few parsecs of the Sun, in which the solar system is embedded (Bruhweiler 



19961) . 



The observation of the stellar continuum below the Lyman edge provides an opportunity 
to directly compare the observed flux distribution of a hot star with the extreme ultraviolet 
(EUV) flux predicted by model stellar atmospheres. Model atmospheres of hot stars are generally 
in good agreement with the observed fluxes longward of 912 A. The lack of spectroscopic data 
in the extreme ultraviolet has, until very recently, left models untested in the EUV. Thus, the 
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observations that the continua of e CMa and (3 CMa observed by EUVE exceed the fluxes predicted 



by the present model atmospheres, by a factor of 5 in the case of /3 CMa (Cassinelli et al. 1996) 



and by more than an order of magnitude in the case of e CMa (Cassinelli et al. 1995, Najarro 



et al. 1996), pose a serious challenge to such computations. 



These stars exhibit an "EUV excess" relative to the predictions of the hydrostatic plane- 
parallel line blanketed LTE models of Kurucz (1992) and the non-LTE models of Pubeny fc Lanz 



(1995) . The non-LTE spherically extended atmosphere models of Najarro et al. (1996), that 
include stellar winds, but neglect metal line blanketing, also do not correctly predict the observed 
EUV flux. The non-LTE atmosphere models of Schaerer &: de Koter (1997), which include stellar 
winds and line blanketing predict the EUV flux distribution of e CMa and (3 CMa, but neglect line 
broadening, which artificially restricts the blanketing effects. When line broadening is included, 
Schaerer &: de Koter find a significant decrease in the predicted EUV flux. These models fail to 
predict the EUV flux presumably because they do not correctly predict either the temperature 
structure or level populations of H and He I in the continuum forming layers. It has also been 
suggested that the temperature structure and level populations are affected by the presence of the 



stellar wind ( [Najarro et al. 1996 , Schaerer & de Koter 1997) and by X-ray heating (Cohen et al 



19961) 



We present here hydrostatic, spherical, non-LTE, metal line blanketed model atmospheres 
and synthetic spectra which closely reproduce both the observed EUV continuum and the flux 
distribution longward of 912 A to 25 /xm of e CMa. In section 2 we describe the calculation of 
our model atmospheres. In section 3 we compare the synthetic spectra to observational data. We 
discuss our results and conclusions in section 4. 



2. Stellar Atmosphere Models 



We have computed stellar atmosphere models with the generalized stellar atmosphere 
computer code PHOENIX (version 8.1). It calculates LTE and non-LTE, line blanketed, spherical 
expanding model atmospheres ( pauschildt et al. 199^ Hauschildt et al. 1996, Hauschildt, 



Baron, fc Ahard 199"7| , [Baron fc Hauschildt 1997| and references therein). Although PHOENIX is 



able to treat stellar winds, we do not consider the effect of a stellar wind here. In addition, we do 
not consider here the effects of rotation. X-ray heating or magnetic fields. PHOENIX can handle 
very large non-LTE model atoms as well as line blanketing by millions of atomic (and molecular) 
lines. This code is designed to be very flexible. It has been used to compute model atmospheres 
and synthetic spectra for novae, supernovae, white, brown, and M dwarfs, and accretion disks in 
Active Galactic Nuclei. The code is very well tested for plane-parallel and static cool stars ( [Allard 
fc Hauschildt 1995[ ). In the following paragraphs we give a brief description of its most salient 
features; more detailed descriptions can be found in the references. 
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We use here an equation of state including up to 26 ionization stages of 39 elements. The large 
number of ionization stages is required to span the wide range of electron temperatures and gas 
pressures encountered in OB-star atmospheres. The models also include all relevant bound-free 



and free-free processes with the cross sections given in the compilations of Mathisen (1984) and 



Verner & Yakovlev (1995). The low densities and complicated radiation fields in the atmospheres 



require that the most important species be treated with self-consistent multi-level NLTE using 



the methods described in [Hauschildt (1993)| , [Hauschildt &: Baron (1995)| , and Hauschildt, Baron, 



& Allard (1997). Line blanketing is a fundamental effect in the computation of OB star model 
atmospheres. Therefore, we include, in addition to the NLTE primary and secondary lines, the 
most important lines dynamically selected from the list of 42 x 10^ lines of Kurucz (1994)| ; for 



details of the selection process and the treatment of these lines see Hauschildt et al. (1995), 
Allard & Hauschildt (1995) and references therein. All models assume solar abundances ( [Anders 
fc Grevesse 19891). 



In order to allow for line scattering in LTE background lines, which is important in low 
density plasmas, we parameterize the albedo for line scattering using a parameter a ( [Hauschildt 



[et al. 1992| ). A fraction, a, of the line absorption coefficient is assumed to contribute to the 
absorption coefficient while the remaining fraction (1-q) is included as part of the line scattering. 
This is done only for lines of species not treated explicitly in NLTE. For our NLTE models we use 
a = 0.95. In our LTE models we use either a = 0.0 for no line scattering, pure LTE, or we use a 
= 0.95 to include line scattering. The effects of LTE line scattering on the temperature structure 
and synthetic spectra are discussed in section 4. 

We include the effects of statistical (random) velocity fields using a Gaussian broadening 
velocity ^. For the models presented here, we use ^ = 2kms~^. Each line is calculated with its 
individual depth dependent intrinsic profile. For strong lines we use Voigt profiles, while for weak 
lines we use Gaussian profiles ( pchweitzer et al. 1996 ). 



The standard outer boundary conditions for our model atmospheres are a continuum optical 
depth of ''"2QQQ^ = 10~^° at 2000 A and an outer gas pressure of Pgas = 10~^ dynes cm~^. The 
inner optical depth boundary is '''^qqq^ ~ 

10^. Models are computed on an optical depth grid of 50 
points. The effects of these boundary conditions and other model parameters on the temperature 
structure and synthetic spectrum are discussed in detail in section 4. 

Our solution of the radiative transfer equation is based on an operator splitting method with 



an exact band-matrix form of the approximate A-operator, see Hauschildt (1992) and Hauschildt 



Storzer, &: Baron (1994) for details. Our numerical solution of the radiative energy equation is 
based on an improved version of the Unsold-Lucy method as described in Allard & Hauschildt 
(1995). The temperature structure is computed from the condition of radiative equilibrium. The 
temperature corrections are based on the absolute error in both the flux and the flux derivative. 
Models are typically converged to better than 1.5% in the absolute flux and the flux derivative at 
every depth. The maximum relative change in the temperature at any depth is not greater than 
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0.8% and better than 0.2% in the EUV continuum forming layers. 

We have computed four different classes of models: LTE models, NLTE models with H and 
He but without metal line blanketing, NLTE models including ions of H, He, C, N, O, and NLTE 
models including ions of H, He, C, N, O, Ca, Mg, and Fe. These classes of models are described 
below. Additional detailed information on the NLTE model atoms can be found in Hauschildt 
et al. (1995), Hauschildt et al. (1996), Hauschildt, Baron & Allard (1997), Hauschildt et al. 
(1997), and references therein. The models discussed in this paper are listed in Table 1. 

1) Typical LTE models (and NLTE models which consider background line blanketing in 
LTE) use lines from ions of H, He, Li, Be, B, C, N, O, F, Ne, Na, Mg, Al, Si, P, S, CI, Ar, 
Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Ga, Kr, Rb, Sr, Y, Zr, Nb, Ba, and La selected 



from Kurucz (1994) , LTE lines are selected for inclusion in the model atmosphere based on the 
strength of the line absorption coefficient at line center, k\c, relative to the absorption coefficient 
of surrounding continuum, Kcont- For most of the models considered here, selected LTE lines 
satisfy Kic/Kcont > 10~^. Under this condition, approximately 2.1 x 10^ lines are considered, of 
which most are treated as Gaussian profiles, while the 9.4 x 10^ strongest lines (Kic/Kcont > 10^) 
are calculated with Voigt profiles. In section 4 we discuss the effects of the number of blanketing 
lines on the model solution. The model is computed on a grid of 11074 wavelength points. The 
LTE models converged using higher resolution grids are discussed in section 4. The LTE models 
are calculated starting from an LTE-grey model and iterated to convergence in 30 iterations, 767 
CPU seconds per iteration, on a HP 9000/770 workstation. 

2) An NLTE pure H-He model with no metal line blanketing uses 101 levels, producing 545 
transitions from the ions H I and He I-II. The models are computed on a grid of 13569 wavelength 
points. The models are started from an LTE-grey model structure and iterated to convergence in 
30 iterations, 782 CPU seconds per iteration, on an HP 9000/770 workstation. 

3) A typical NLTE H-He-CNO model includes 1751 distinct levels, producing 11975 transitions 
from the ions H I, He LII, C I-IV, N I- VI, and O LVI in detailed NLTE in addition to the 
background LTE line blanketing discussed above. The models are computed on a grid of 63486 
wavelength points. The NLTE models are calculated starting from a converged LTE line blanketed 
model and iterated to convergence in 30 iterations, 7707 CPU seconds per iteration, on an HP 
9000/770 workstation. 

4) More complete NLTE models employ 3035 levels, with 37151 transitions from the ions H 
I, He LII, C LIV, N I- VI, O I- VI, Ca II, Mg II, and Fe II-III in NLTE in addition to the LTE 
line blanketing. These models are computed on a wavelength grid of 183591 wavelength points. 
The models are calculated starting from a NLTE H-He-CNO class model and converged in 10 
iterations, 8900 CPU seconds per iteration, running the parallel version of the PHOENIX code on 
an IBM SP2 computer at the Cornell Theory Center. 



- 6 - 



3. Comparison of Models with Data 



We compare here synthetic spectra generated from our model atmospheres with the observed 
spectrophotometry. The EUVE data are from Cassinelh et al. (1995). The archival UV data 
are from the S2/68 Ultraviolet Sky Survey Telescope (UVSST) on TD-1 (|Jamar et al. 1976D , 
the Orbiting Astronomical Observatory 2 {OAO-2 ) ( pode &: Meade 1979D , the International 
Ultraviolet Explorer (lUE ), and the Goddard High Resolution Spectrograph (GHRS ). Optical 
spectrophotometric data are from Johnson fc Mitchell (1975) and Burnashev (1985) . The TD-1 
and OAO-2 data were obtained from the on-line archival catalogs of the Astronomical Data Center 
(ADC) at Goddard Space Flight Center (GSFC). The lUE data were obtained from the National 
Space Science Data Center at GSFC. With the exception of the long wavelength primary (LWP) 
camera high-dispersion spectra, the lUE data are all processed with the NEWSIPS software ( |De| 
La Pena et al. 1994] ). The GHRS data ( [Heap et al. 1995 ) were obtained from the GHRS data 
reduction facility at GSFC and reduced using the CALHRS software in IDL maintained by the 
GHRS Science Team. Infrared J and K band fluxes are from Johnson et al. (1966) and the IRAS 
color-corrected fluxes at 12 /xm and 25 //m are from Cassinelli et al. (1995). A synthetic spectrum 
from model Nl is compared with the observed spectrum of e CMa from 350 A to 25 fim in Figure 
1. The stellar parameters of model Nl are listed in Table 1. 



3.1. EUV Data 



Synthetic spectra from models Nl and N2 are compared with the observed EUV spectrum 
of e CMa in Figure 2. A plane-parallel LTE ATLAS9 model {Tes = 21000 K, log 5 = 3.2, ^ 
= 2 kms"^) (^ urucz 1992| ) is shown for comparison. The effective temperature of models Nl 
and N2, Teg = 21750 K, is at the limit of the range established by Code et al. (1976): Tefj = 
20990ib760K. The model surface gravity, logg = 3.2, and angular diameter, 9^ = 0.77 mas, are 
within the established limits: logg = 3.20±0.15, 9^ = 0.80±0.05 mas (Cassinelli et al. 1995). 
Comparison of the position of e CMa on the HR diagram with evolutionary tracks suggests the 
surface gravity is closer to logg = 3.5 (see section 4). The synthetic spectra are scaled to match 
the observed optical and near-IR spectrophotometry by adopting an angular diameter of 0.77 mas 
The ATLAS9 model is scaled in the same way. We emphasize here that the data have not been 
normalized - all comparisons are made based on absolute fluxes. 



Fig. 1. — The continuous energy distribution of e CMa from 350 A to 25 ^m compared with the 
21750 K NLTE model (Nl) and the ATLAS9 21000 K model. The EUVE data are corrected for 
the interstellar neutral hydrogen column assuming a neutral hydrogen column density of N(H I) = 
5 X 10^^ cm-2. 
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The EUVE absolute calibration uncertainty, which awaits a complete in-orbit analysis, is 
estimated to be 25%, while the relative error is better than 10% ( |Vallerga fc Welsh 1996|) . We 



correct the EUVE data for transmission through the neutral interstellar column based on cross 
sections of Rumph, Bowyer, & Vennes (1994) and Morrison & MacCammon (1983), using the 



on-line ISM Transmission Tool available at the Center for Extreme Ultraviolet Astronomy Web 
Page. We adopt a neutral hydrogen interstellar column density of N(H I) = 5 x lO^^cm"^, the 
upper limit established by Gry et al. (1995). 

Our synthetic spectra of models Nl and N2 very closely reproduce the H I Lyman continuum 
observed by EUVE. Model N2, which considers ions of H, He, and CNO in NLTE, provides the 
best fit to the data below the 504 A He I edge. However, the more physically realistic model Nl, 
which includes Mg II, Ca II, and Fe II-III in NLTE, predicts less flux in the He I continuum. At 
400 A the Nl synthetic spectrum is a factor of 3 below the observed flux. 

A synthetic spectrum with the identical temperature structure of model Nl, but without 
metal line blanketing, shows that the shape and flux level of the H I and He I Lyman continua 
are not severely distorted by metal line blanketing (see Figure 3). At 1 A resolution, comparable 
to that of the EUVE spectrum, the level of the apparent continuum is lowered by at most 20% 
between 500 A and 700 A. 

The temperature structures and EUV synthetic spectra of models Nl, N2, and N7 are 
compared with an LTE model, LI, in Figure 4. The H I Lyman continuum in both the NLTE 
and LTE models forms at column masses greater than 10"^ g cm^^ on the plateau portion 
of the temperature structure common to these four models at approximately 17200 K. The 
overpopulation of the H I ground-state in the NLTE models relative to the LTE model causes the 
NLTE models to have a lower EUV flux immediately below 912 A due to the larger cross section 
for H I ionization from the ground-state. Far from the 912 A edge, just above the 504 A He I edge, 
where the H I ionization cross section approximately an order of magnitude smaller, the NLTE 
and LTE fluxes are in better agreement. NLTE appears not affect the temperature structure in 
the H I Lyman continuum forming layers. The H I departure coefficients, however, regulate the 
flux in the Lyman continuum. The NLTE models are significantly cooler at column masses less 
than 10~^ g cm~^ compared with the LTE model such that the He I edge, which is absent in the 
LTE model spectrum, is quite strong in the NLTE models. The larger relative fraction of He I 
and the overpopulation of the He I ground-state in the NLTE models results in a significantly 
lower EUV flux below 504 A compared with the LTE model. An ATLAS9 model spectrum with 
parameters Tefj = 21750 K and \ogg = 3.2 is shown in Figure 4 for comparison. Both the LTE 



Fig. 2.— The EUVEIM spectrum of e CMa from 350 A to 700 A compared with the 21750 K NLTE 
models (Nl, N2, and N7) and a ATLAS9 21000 K model. The EUVE data are corrected for the 
neutral interstellar column adopting a neutral hydrogen column density of N(H I) = 5 x 10^'^ cm~^. 
The PHOENIX model synthetic spectra have a resolution of 1 A. 
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and NLTE spherical PHOENIX models predict a higher EUV fluxes than the ATLAS model. This is 
primarily due to the significant differences in outer temperature structures between spherical and 
plane-parallel models. We discuss these differences in detail in section 4. 



3.2. Ultraviolet Data 



The extremely low column density, N(H I) < 5 x 10^^ cm~^ ( |Gry et al. 199^ ), toward e 
CMa produces negligible dust extinction in the ultraviolet and visible portions of the spectrum. 
Consequently, no such corrections were applied to these data. 

The pre-COSTAR GHRS large science aperture (LSA) ECH-B data provide a measurement of 
the absolute UV continuum of e CMa at four 15 A wide orders between 2350 A and 2850 A. Figure 
5 shows the comparison between a synthetic spectrum from model Nl and the UV and optical 
spectrophotometry. We plot the data and synthetic spectrum in Figure 5 as A versus X'^Fx to 
accentuate small differences in the slope and absolute levels of these spectral energy distributions 
in the ultraviolet region. The absolute UV flux measurements from the GHRS ECH-B spectra 
are in good agreement with model Nl. The GHRS spectra and model Nl show ~30% more flux 
than measured by either lUE, OAO-2, or TD-1. The difference exceeds the published errors in 
the absolute calibration of these data. The published uncertainty for TD-1 is ±20% and the 
uncertainties for the OAO-2 data are 30% below 1300 A, 10% between 1300 A and 1800 A, and 
longward of 1800 A ( peeckmans 1977 ). The established uncertainty in the absolute calibration 



for lUEis ±10% (iBohlin 1986|) 



In light of the agreement between model Nl and the GHRS, optical, and IR fluxes for e CMa, 
the lUE, TD-1, and 0^0-^ data appear to have a wavelength independent systematic calibration 
error of ~30%. The OAO-2 data (|Code fc Meade 197g| ) do not include correction of Bohlin k 



Holm (1984) based on the choice of r] UMa (HR 5191, HD 120315; Tes = 17000 K, logg = 4.25) 
as a fundamental UV standard star. The UV spectral energy distribution of r] UMa, based on 
multiple lUE observations, when compared with the OAO-2 data shows that the OAO-2 flux is 
systematically larger between 1280 A and 2040 A, peaking at ~ 30% larger around 1550 A. Bohlin 



fc Holm (1984) correct for this OAO-2 Hux excess. Applying this same correction to the OAO-2 



data for e CMa improves the agreement between these data and the lUE and TD-1 data sets. 

Only a single scan of e CMa was made with Spectrometer 1 (range : 1850 A - 3600 A) by 
OAO-2 ( Pode &: Meade 1979 ). Two scans were obtained with Spectrometer 2 (range: 1160 A 



1850 A). Spectra of stars for which we find good agreement between our model absolute flux and 



Fig. 3. — Model N2 in the EUV region with and without metal-line blanketing. The synthetic 
spectra have a resolution of 1 A. 
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the OAO-2 absolute flux, such as a Vu', consist of six scans in Spectrometer 1 and three scans in 
Spectrometer 2. Since the OAO-2 absolute flux of e CMa longward of 1850 A is based on a single 
scan it may be subject to greater error. In addition, the (9^ integrated absolute flux for e CMa 
from 1810 A to 3300 A as listed by Code et al. (1976) is 20% greater than the value we compute 
from the Code &: Meade (1979) data set. For the integrated flux of a Vir in the same wavelength 
region we flnd no such disagreement between the two data sets. So, the Code et al. fluxes which 
are used to compute the effective temperature are in better agreement with the GHRS absolute 
fluxes than the OAO-2 data shown in Figure 5. The absolute flux of the spectrophotometric 
standard jJL Col (HR 1996, HD 38666) at 2740 A as measured by GHRS (3.1 xlO'^^ erg cm'^ s'^ 
A~^; Heap et al. 1995 and references therein) and OAO-2 (2.83xl0~^° erg cm" -2 A-i; Code 
& Meade 1979) differ by just less than 10%. Like e CMa, a single scan of Col was obtained by 
Spectrometer 1. 

A discrepancy in the relative flux distribution was found between model atmospheres and the 
0^10-^ data for [5 CMa ( [Cassinelh et al. 1996| ). Applying the correction of [Bohlin k Holm (1984J 



to the OAO-2 data of fi CMa appears to resolve the discrepancy in the relative flux distribution 
for this star. In this case, however, neither GHRS nor lUE large aperture data longward of 2000 A 
are available with which to check the absolute calibration. As with e CMa, the only a single scan 
of (3 CMa was obtained by Spectrometer 1. We find the Code et al. (1976) absolute integrated 
flux between 1850 A and 3300A to be 10% larger than the same measurement from the Code & 
Meade (1979) data set. The absolute flux from the lUE large aperture SWP data for /3 CMa is in 
agreement with our model atmosphere for this star (Aufdenberg et al. 1997) . 

Here we investigate possible reasons for the discrepancy in the absolute flux level between the 
model spectra and the lUE observations. Low-dispersion, large aperture, lUE observations of a 
star as bright as e CMa (V = 1.5) required a rapid trailed exposure. Only two such exposures of 
e CMa exist, SWP54338 and LWP30419, with exposure times of 0.184 s and 0.202 s respectively. 
Rapid trails with exposure times less than a second, are generally not as reliable as slower trailed 



spectra (Imhoff 1997). For example, the star may not have crossed the lUE large aperture along 
the nominal path, thus affecting the overall flux levels. We inspected the low dispersion re-sampled 
lUE images (SILO) images of e CMa and compared them to a rapid trailed SILO image of a Vir 
(SWP39248, 0.179 s), for which there is very good agreement (within 5%) between the measured 
absolute flux and our model atmosphere (Aufdenberg et al., 1997). The e CMa SWP and LWP 
spectra are noticeably narrower than the width of the a Vir spectrum perpendicular to the 
dispersion axis. This may account for the 30% discrepancy in the lUE data. An additional source 
of error in the absolute calibration is the extrapolation of the SWP and LWP camera sensitivity 
degradation functions outside beyond valid date range. Correcting for this increases the absolute 



Fig. 4. — Temperature structures (left) and synthetic spectra in the EUV (right) for models Nl, 
N2 and N7. A spherical LTE model (LI) and an ATLAS9 spectrum are included for comparison. 
The PHOENIX synthetic spectra have a resolution of 1 A. 
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flux by ~4% (Imhoff 1997). 



3.2.1. lUE High Dispersion Data 



The lUE low dispersion spectra were used by NEWSIPS to establish the absolute calibration 
applied to the large aperture high resolution spectra. This explains why the lUE echelle spectra 
also disagree with the model and the GHRS fluxes. Despite the absolute flux discrepancy between 
the lUE data and the model, there is good agreement with the observed relative flux distribution. 
A comparison of model Nl and lUE high dispersion spectra SWP54337 and LWP30183 is shown 
in Figure 6. These high dispersion images were chosen for their well exposed and non-saturated 
continua. Non-saturated regions of the other high dispersion spectra we have examined show the 
same absolute flux level. 



3.2.2. GHRS Data 



The model Nl flux agrees with the GHRS continuum flux to within ±2%. In addition to 
the agreement between the model and the absolute flux level of the four GHRS ECH-B orders, 
high resolution synthetic spectra of model Nl are in reasonable agreement with the observed line 
strengths and profiles (see Fiugres 7-9). The high resolution model spectra have been scaled by a 
small factor (< it 2%) relative to the overall fit for the continuous energy distribution to match 
the observed continuum level. The widths of the observed line profiles are slightly increased as a 
result of the spherical aberration in these pre- COSTAR observations Peap et al. 1995| ). The 
synthetic spectra were rotationally broadened to 35 kms~^ ( Passinelli et al. 1995 ) and corrected 
for the heliocentric radial velocity of e CMa, +27.4 kms~^( |Gry et al. 1995 ). No attempt was 
made to improve the match between the observed and synthetic spectra by changing the chemical 
abundances or other model parameters. 



3.3. Infrared Data 



Fig. 5. — Comparison of UV and optical spectrophotometry to model Nl. With the exception of 
the GHRS LSA fluxes, the observed UV spectra lie below the model spectrum by as much as 30%. 
The correction of Bohlin & Holm (1984) has not been applied to the OAO-2 data shown here. The 
synthetic spectrum has a resolution of 5 A. 
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Comparison of model Nl mid-IR fluxes with the color-corrected IRAS fluxes at 12 /xm and 25 

//m from Cassinelli et al. (1995) indicates the presence of a slight mid-IR excess. The IRAS fluxes 
are 7% larger at 12 fim and 9% larger at 25 fim compared with model Nl. Model Nl fluxes are 
within 1.1 a and 1.2 a of the mean IRAS fluxes at 12 //m and 25 fim, respectively. The model 
predicts 9% and 10% more flux at 12 and 25 fim respectively, compared with the ATLAS9 
Teff = 21000 K, log 5 = 3.2, 6*^ = 0.80 mas model. 



4. Discussion and Conclusions 

4.1. Stellar Parameters and the Local ISM 

We flnd that the observed flux distribution between 350 A and 25 fim can be reasonably well 
represented by our spherical, hydrostatic, non-LTE, line blanketed model atmosphere by adopting 
model parameters within the range of uncertainty of the established stellar parameters. Adopting 
the HIPPARCOS distance of 132^g^ pc and the best fit model parameters, 6d=0.77 and logg 
= 3.2, yield a radius of 11.3ib 1.1 Rq and a mass of 7.01q'9 Mq for e CMa. The uncertainties 
in radius and mass reflect only the uncertainty in the HIPPARCOS parallax. Our value for the 
radius of e CMa is ^ 30% smaller than that quoted by Cassinelli et al. (1995) primarily due to 
the HIPPARCOS distance, which puts e CMa 30% closer than previously believed. This smaller 
radius reduces the mass by 54% for the same surface gravity. The surface gravity of logg' = 3.2 
appears to be too low by approximately 0.3 dcx. Based on a comparison with evolutionary tracks 
(sec Figure 10), e CMa has a mass of nearly 12Af0, however, the measured radius with a gravity 
of logg = 3.2ib0.15 yields a mass of 7.4ib3.0MQ. An atmosphere model (model N7) computed 
with logg = 3.5, but with otherwise identical stellar parameters as model Nl, does not show a 
significantly diflFerent EUV flux distribution (see Figure 2). A possible resolution to this mass 
discrepancy may result from a reanalysis of the Balmcr line flts used to establish the surface 
gravity. Apparently the only published high dispersion Balmer proflles are photographic data from 
Buscombe (1970). We leave this analysis for future work. 

The predicted number of hydrogen and helium ionizing photons at the stellar surface from 
model Nl are 1.59 x 10^^ photons cm"^ s"^ (0 A - 912 A) and 1.13 x 10^^ photons cm^^ s"^ 
( A - 504 A) respectively. The best fitting model EUV fiux (model N2) at the stellar surface 
is listed in Table 2. Using the best fit angular diameter, 0.77 mas, the calculated unattenuated 
hydrogen ionizing flux (0 A - 912 A) and neutral hydrogen photoionization rate at the Local Cloud 

Fig. 6. — lUE high-dispersion data compared with model Nl. The data have been scaled by a 
factor of 1.36 to match the model flux. The lUE data have been binned to 1 A to match the 
synthetic spectrum resolution. 
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surface are 2290 photons cm~^ s~^ and 1.06 x 10~^^ respectively. These values are consistent 

with those of Vallerga &: Welsh (1996) who compute values more than twice as large but adopt a 
correspondingly larger neutral hydrogen column density of 9 x 10^^ cm~^ and a distance of 188 pc. 



4.2. Effects of Sphericity and Line Blanketing on the EUV Continuum 

In this section we explore the sensitivity of our model atmospheres to the treatment of 
the problem with spherical or plane-parallel geometry, the degree of line blanketing, and other 
model parameters. In Figure 4 we show synthetic spectra from our LTE and NLTE spherical 
line blanketed models compared with a LTE plane-parallel line blanketed model from Kurucz 
with same effective temperature and surface gravity. We find that the differences in the predicted 
EUV flux between our models and that of Kurucz are primarily due to the differences in the 
temperature structures between line blanketed spherical models and line blanketed plane-parallel 
models. To explore these differences, we computed a grid of models, all of which have an effective 
temperature of Tcfj = 20990 K, a surface gravity of log 5 = 3.2, and a turbulent velocity of ^ = 2 
km s~^. These models differ in the geometry of the calculation (spherical or plane-parallel), the 
use and degree of LTE or NLTE line blanketing, the presence or absence of LTE line scattering, 
the value of the outer boundary pressure, the value of the outer boundary continuum optical 
depth, and the resolution of the wavelength grid in the model calculation. These models are listed 
in Table 1. 

We find that the combination of parameters in our models which has the most significant 
effect on the temperature structure, and therefore the EUV flux level, is the inclusion of both metal 
line blanketing and a spherical geometry. The temperature structure, in particular temperatures 
in the region responsible for the EUV continuum, is only slightly affected, if at all, by different 
hydrostatic outer boundary conditions, the degree of line blanketing, the presence of LTE line 
scattering, or the resolution of the wavelength grid. We now discuss briefly each of these test 
models. 



4-2.1. Sphericity and Line Blanketing 

In Figure 11 we compare the temperature structures of six LTE and NLTE models with 
and without line blanketing. The ATLAS9 temperature structure is shown for comparison. The 

Fig. 7. — GHRS ECH-B observation Z14Z050MT compared to high resolution synthetic spectrum 
of model N2 
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Fig. 8. — GHRS ECH-B observation Z14Z050XT compared to high resolution synthetic spectrum 
of model N2. 



Fig. 9. — GHRS ECH-B observation Z14Z051BT compared to high resolution synthetic spectrum 
of model N2. 



Fig. 10. — Error box for e CMa based on the observational uncertainties plotted along with 
evolutionary tracks (Schaller et al. 1992) of solar metallicity on a theoretical HR diagram. The 
error bars reflect the combined 1 a errors in the measured total flux and the measured angular 
diameter (Code et al. 1976), and the standard error in the HIPPARCOS parallax. The diamond 
symbol represents the position of e CMa from the best fit values based on the comparison of the 
synthetic spectra to the spectrophotometry adopting the mean HIPPARCOS distance. 



Fig. 11. — Temperature structures for a variety of spherical and plane-parallel models listed in 
Table 1. An ATLAS9 structure is included for comparison. Spherical line blanketed models show 
significantly warmer temperature structures compared with plane-parallel line blanketed models. 
NLTE models with no metal line blanketing show identical structures for the two geometrical 
treatments. 



Fig. 12. — Synthetic spectra in the EUV for a variety of spherical and plane-parallel models listed in 
Table 1. An ATLAS9 spectrum is included for comparison. Spherical line blanketed models predict 
significantly larger EUV fiuxes compared with plane-parallel line blanketed models. NLTE line 
blanketed models predict similar or lower EUV fluxes than LTE models due to the overpopulation 
H I and He I ground states in the NLTE treatment. NLTE models with no metal line blanketing 
show identical EUV spectra for the two geometrical treatments. 



-14- 



temperature structures of the line blanketed models fall into two groups. The spherical models 
show a plateau structure near 16600 K between column masses 1 g cm~^ and 10~^ g cm~^, while 
the plane-parallel models show a plateau structure near 13300 K between column masses 10^^ g 
cm~^ and 10"*^ g cm~^. The NLTE and LTE models have nearly the same temperature structure 
in these plateau regions and differ substantially only beyond the plateau at lower mass columns. 
The differences in the temperature structures between the two geometries have significant effects 
on the flux levels of the synthetic EUV continua as shown in Figure 12. For example, the EUV 
flux at 750 A is a factor of 6 higher for the spherical LTE line blanketed model (L2) than the 
plane-parallel LTE line blanketed model (L3). 

The H I Lyman continuum of the LTE spherical model (L2) forms at 16500 K, while for the 
LTE plane-parallel model (L3) the Lyman continuum forms at approximately 14500 K. The Lyman 
continua of the NLTE line blanketed models form at the same column masses in the atmosphere 
as the LTE models, but, the overpopulation of the ground states in the NLTE models further 
affects the flux levels of the H I and He I Lyman continua. Our LTE, line blanketed, plane-parallel 
synthetic spectrum predicts the same flux in the H I Lyman continuum as the ATLAS9 spectrum, 
but these models predict different flux levels for the He I continuum below 504 A. This is due 
to the different temperature structures predicted by these models for column masses less than 
10~^ g cm~^. We do not yet understand the reasons for the differences between these models in 
this region. Differences in the two codes may exist in the treatment of metal line blanketing. In 
addition, important differences may exist between the two codes for the equation-of-state, which 
is computed by PHOENIX in a completely independent way from that of ATLAS, e.g., by using 
different energy level data. 

In contrast to the large differences in the temperature structures and EUV continua for the 
line blanketed models with respect to a spherical versus a plane-parallel geometry, the NLTE 
H-He only models (N5 and N6), with no metal line blanketing, show nearly identical temperature 
structures and EUV spectra for both geometries. Therefore, the significant differences in the 
model EUV spectra of early B stars for a particular treatment of the geometry is the result of the 
inclusion of metal line blanketing. The Lyman continuum of the unblanketed pure H-He NLTE 
model forms at 13500 K and the flux level at 750 A is more than a factor of 15 lower than that 
predicted by the spherical NLTE model (N3) with line blanketing. 

We suggest that this can be explained physically as follows. Including metal line blanketing in 
a model atmosphere produces the well known effects of backwarming and surface cooling relative 
to the temperature structure in an unblanketed atmosphere. These effects are seen by comparing 
the temperature structures of plane-parallel model, L3 (metal line blanketed) with model N6 (H 
and He opacity only) in Figure 11. The surface layers, mass depths < 10~^ g cm~^, of the L3 
structure are cooler relative to N6, while the deeper layers, depths > 10~^ g cm~^, are warmer 
relative to N6 until depths greater than 1 g cm~^ are reached. In the layers of the atmosphere 
sensitive to blanketing, the metal lines produce an opacity, primarily in the ultraviolet, that heats 
the layers in radiative equilibrium immediately below and cools the layers immediately above. 



-15- 



producing a steeper temperature gradient relative to the model lacking this metal opacity. 

These backwarming and surface cooling effects of line blanketing are seen in the spherical 
models, but with a difference. Looking at the temperature structures of models L2 (spherical) 
and L3 (PP), these models have basically the same slope from 10~^ g cm~^ up to their respective 
plateaus. The spherical and PP structure shapes are very similar from the plateaus outward. 
However, the PP model temperature structure flattens out at 10~^g cm~^, while the spherical 
model temperature slope is sustained into deeper layers (to 10~^ g cm~^) before flattening out 
to its plateau. The spherical temperature structure maintains the same slope, from 10~^ g cm~^ 
inward to the plateau, as the PP case, but to deeper layers. In other words, the temperature 
gradient outward from the plateau begins at greater depth in the spherical model compared with 
the PP model. Effectively, in the spherical case, the surface cooling begins "earlier", moving 
outward in the structure. Therefore in the spherical case, the resultant backwarmed layers (depths 
> 10~^ g cm^^) which form the plateau, are at a higher temperature beneath the cooling layers 
(depths < 10~^ g cm~^) relative to the plateau in the PP model, where cooling begins further 
outward. Since the surface cooling begins at a greater depth in the spherical case, the backwarming 
is shifted to deeper layers relative to the PP case. These deeper layers are intrinsically warmer, so 
the temperature of the spherical model backwarmed plateau sits at a higher temperature than the 
PP model plateau. 

One important factor which shifts the surface cooling to greater depths in the spherical 
models is the reduced optical depth seen by photons in a spherical geometry relative to that in a 
semi-infinite plane-parallel geometry at a given physical depth. Immediately below the cooling 
layers, the enhanced photon escape probability increases the populations in the lower levels of the 
metal transitions, thereby increasing the local opacity. In response, the layers immediately below 
must heat up in order to transport the same flux in the presence of the increased line opacity 
and reduced frequency bandwidth. As a result, in the spherical models, backwarming begins at a 
greater depth than in the PP models. 

The result that the H-He only models (N5 and N6) arc not sensitive to the geometry is likely 
because the H-He line opacity, which is nearly absent in the ultraviolet, has a much weaker effect 
on the temperature structure relative to the metal lines. The enhanced escape probability for 
photons in the spherical geometry of a H-He atmosphere has a much weaker effect on the cooling 
rates and therefore the temperature structure. 

Our results can be explained on the basis of the geometrical effects in combination with the 
metal line blanketing in LTE. NLTE processes, which are necessary for the detailed modeling of 
the observed spectra (for instance, the presence of the He I 504 A edge), clearly have effects on the 
temperature structure, but these are less significant that those effects due to the combination of a 
spherical geometry and LTE metal-line blanketing alone. 
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4-2.2. Number of LTE Blanketing Lines 

We have tested the sensitivity of the temperature structure to the number of blanketing atomic 
lines. For the typical LTE model the number of lines used in the calculation is approximately 
2 X 10^. The temperature structure, however, is mainly determined by approximately the 50000 
strongest lines. Figure 13 compares the temperature structures and synthetic spectra in the EUV 
for a grid of six spherical models (L2, L11-L15) each with a different number of atomic blanketing 
lines from 2095181 lines to only the strongest 217 lines. 

The level of the plateau in the temperature structure, between column masses 10^^ g cm^^ 
and 10~^ g cm~^, is 600 K lower for 217 lines compared with 2095181 lines. The level of the 
EUV continuum at 750 A is 30% lower for 217 lines compared with 2095181 lines. In contrast, a 
line blanketed plane-parallel model (L3) is 2800 K cooler than an otherwise similar line blanketed 
spherical model and their predicted EUV fluxes differ at 750 A by more than a factor of 6. We 
conclude the EUV continuum flux in LTE is regulated by the several thousand strongest blanketing 
lines and that there is little affect on this EUV continuum from the inclusion of millions of weak 
lines used in the model. 



4-2.3. LTE Line Scattering 

In order to test the effect of LTE line scattering on the temperature structure of LTE line 
blanketed models, we have computed both spherical and plane-parallel models (L2-L5) with 
values for the depth independent line scattering albedo of a = 95% and a = 0%. Neglecting line 
scattering corresponds to pure LTE. The temperature structures and EUV synthetic spectra are 
compared in Figure 14. Including LTE line scattering reduces the temperature down to column 
masses of 1 g cm~^ and consequently also reduces the flux level of the EUV continuum. The line 
scattering in our LTE models changes the level of the EUV continuum by factors of ~1.5 in the 
EUV continuum at 750 A and ^2.5 at 450 A. These effects are less significant than those due to a 
plane-parallel geometry. 



4-2-4- Outer Boundary Conditions 

Fig. 13. — Temperature structures (left) and synthetic spectra in the EUV (right) for models which 
consider different numbers blanketing lines in LTE. The temperature structure and EUV spectrum 
are sensitive to approximately the strongest 10,000 lines and a largely unaffected by millions of 
weaker metal lines. 
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Our standard value for the outer boundary gas pressure, Pgas = 10~^, is one which is 
traditional (Kurucz 1970). To test the effect of this assumption wc have converged models (L2, 
L6, L7) for outer gas pressures of 10~^, 10~^, and 10~^ dynes cm'^. These models produce nearly 
identical temperature structures and identical synthetic spectra (see Figure 15). We have also 
converged models (L9 and LIO) for outer continuum optical depth of 10~^ and 10~^°. Adjusting 
the outer boundary optical depth from 10""^'^ to 10~^, similar to the ATLAS9 model, appears to 
have no effect on the temperature structure or synthetic spectra for a given geometry, spherical or 
plane-parallel (see Figure 16). Realistic outer boundary conditions for e CMa are probably not 
hydrostatic. We leave the investigation of the effects of more realistic boundary conditions, such 
as a stellar wind, for future work. 



4-2.5. Wavelength Grid Resolution 

Wc have computed a small grid of spherical LTE line blanketed models to test the sensitivity 
of the converged temperature structure to the number of wavelength grid points employed in 
the model calculation. Figure 17 shows the temperature structures and EUV synthetic spectra 
for three models which were computed on wavelength grids of 11074, 22148, and 44296 points 
respectively. The 11074 grid points model has a slightly warmer temperature structure of column 
masses less than 10~^ g cm~^ and a more linear slope near 10~^ g cm~^. At greater depths the 
temperature structures of the three models are nearly identical and the synthetic EUV spectra of 
the three models arc also nearly identical. The NLTE models that were computed for comparison 
with observational data use more than 60000 wavelength points. 



4.3. Summary 

The continuous energy distribution of e CMa from the EUV to the IR is reasonably well 
fit by a hydrostatic model atmosphere that neglects a stellar wind, magnetic fields, and X-ray 
heating. We find the fiux level of the helium continuum to be sensitive to the number of ions and 
transitions treated in NLTE. Our most realistic spherical, NLTE, hydrostatic model is unable 
to reproduce the slope of the observed spectrum below 504 A. The slope may be affected by a 

Fig. 14. — Temperature structures (left) and synthetic spectra in the EUV (right) for models which 
have either a large value, (95%), for the LTE line scattering albedo, or no LTE line scattering. The 
temperature structure is generally slightly cooler with the inclusion of line scattering and the EUV 
continuum flux is slightly weaker. 
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stellar wind and non-thermal heating processes which may alter the ionization structure of the 
atmosphere and the population of the He I ground state. Our model fit, based on absolute fluxes 
together with the HIPPARCHOS distance, yields fundamental parameters for e CMa. 

Our Tefj = 21000 K, log 5 = 3.2, spherical NLTE model predicts more than twice as many 
hydrogen ionizing photons and over 200 times more neutral helium ionizing photons compared 
to a standard hydrostatic plane-parallel LTE model with the same stellar parameters. This is 
important since the EUV flux from e CMa may not be peculiar, and is well represented by our 
spherical hydrostatic model atmosphere. It touches directly on the contribution of early B stars 
to the diffuse UV radiation field and the photoionzation of the interstellar medium. 

We find significant differences in the strength of the predicted EUV flux between our line 
blanketed spherical and line blanketed plane-parallel models. A more realistic model treatment of 
early B giants with a spherical geometry and NLTE metal line blanketing results in the prediction 
of significantly larger EUV fluxes compared with plane-parallel models. This result appears 
to explain a large part of the reported discrepancy between the observed EUV flux of early B 
giant e CMa and the EUV flux predicted by plane-parallel LTE and NLTE line blanketed model 
atmospheres. 
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significantly improved this paper, D.H. Cohen for providing us with the reduced EUVE data and 
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Fig. 15. — Temperature structures (left) and synthetic spectra in the EUV (right) for models with 
different values of the outer boundary gas pressure (in dynes cm~^). The temperature structure 
and EUV flux are little affected by the value of the outer pressure. 
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Fig. 16. — Temperature structures (left) and synthetic spectra in the EUV (right) for models with 
different values of the outer boundary continuum optical depth. The temperature structure and 
EUV flux are not affected by the value of the boundary optical depth, however, they are significantly 
affected by the geometric treatment of the model atmosphere. 



Fig. 17. — Temperature structures (left) and synthetic spectra in the EUV (right) for models 
with different resolutions of the wavelength grid used in the model computation. The temperature 
structure and the EUV flux are not significantly affected by the resolution of the wavelength grid. 
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Table 1: Models 



Modei 


Type 


J-cS (J^j 


log (5) 


^ ( km s ^ ) 


vjeometry 


LIE hcattermg 


Comment 


ATI 

IN i 


ATT 1^T7 




6.2 


Z.O 


bpnerical 


yes 


rl-xie-UiN U-Ua-iVig-r'e 


JNz 


INLlli/ 


zi75(J 


6.2 


Z.O 


bpnerical 


yes 


H-He-UiNU 




ATT T'TT' 


zuyyu 


o.z 


Z.O 


Spherical 


yes 


n-ile-UJNU 


AT /I 

JN4 


ATT T^T^ 


zuyyu 




6.2 


Z.O 


Fiane-r'araiiei 


yes 


n-ile-UlNU 


ATC 

JN5 


ATT T^TT^ 


zuyyu 


6.2 


2.U 


Dpnerical 


no 


n-ile only 


ATt; 

JNd 


ATT T^T? 


zuyyu 


6.2 


Z.O 


T>i „ T ^ „ n „i 

Fiane-Farailei 


no 


1 J 11 1 

H-He only 


ATT 

N7 


ATT 'T~1T7^ 


21750 


3.5 


2.0 


hpnerical 


yes 


TT TT /^AT/^ TV T T7< 

H-He-CJN (J-Ca-Mg-l'e 


LI 


LTE 


21750 


3.2 


2.0 


Spherical 


yes 




L2 


LTE 


20990 


3.2 


2.0 


Spherical 


yes 




L3 


LTE 


20990 


3.2 


2.0 


Plane-Parallel 


yes 




L4 


LTE 


20990 


3.2 


2.0 


Spherical 


no 




L5 


LTE 


20990 


3.2 


2.0 


Plane-Parallel 


no 




L6 


LTE 


20990 


3.2 


2.0 


Spherical 


yes 


P — 10~^ 


L7 


LTE 


20990 


3.2 


2.0 


Spherical 


yes 


Pga. = 10-6 


L8 


LTE 


20990 


3.2 


2.0 


Spherical 


yes 


Tmin = 10-6 


L9 


LTE 


20990 


3.2 


2.0 


Plane-Parallel 


yes 


Tmin = 10-6 


LIO 


LTE 


20990 


3.2 


2.0 


Spherical 


yes 


l^lcl '^cont > 1 


Lll 


LTE 


20990 


3.2 


2.0 


Spherical 


yes 


Klc/Kcont > 10 


L12 


LTE 


20990 


3.2 


2.0 


Spherical 


yes 


Klc/^cont > 10^ 


L13 


LTE 


20990 


3.2 


2.0 


Spherical 


yes 


«;ic/'^cont > 10^ 


L14 


LTE 


20990 


3.2 


2.0 


Spherical 


yes 


Klc/^cont > 10^ 


L15 


LTE 


20990 


3.2 


2.0 


Spherical 


yes 


22148 wavelength points 


L16 


LTE 


20990 


3.2 


2.0 


Spherical 


yes 


44296 wavelength points 
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Table 2: Model N2 EUV Surface Flux 



A (A) 


Log Fa^ a (A) 


T,np- F\^ 


2f)n 


-O.Oli 


oyu 


7 376 






DUO 








DZU 


7 557 


245 


U.OOO 


DoO 


7 616 






DOU 


7 584 


275 


1 7Q^ 


ODO 


7 6Q4 


290 


9 "^98 


OoU 


7.841 


305 


9 8n^ 


Quo 


7.901 


320 


079 


71 n 

1 lU 


7.847 


335 




79 1^ 
1 ZO 


7.979 


350 




7/1 n 


8 054 


3fi5 


4.261 


755 


8 0Q7 


ooO 


4.535 


770 


8.13/ 


395 


4.786 


785 


8.175 


410 


5.040 


800 


8.212 


425 


5.225 


815 


8.245 


440 


5.456 


830 


8.280 


455 


5.654 


845 


8.305 


470 


5.827 


860 


8.334 


485 


5.978 


875 


8.362 


500 


6.138 


890 


8.223 


515 


6.726 


905 


8.409 


530 


6.807 


920 


8.990 


545 


7.140 


935 


8.830 


560 


7.115 


950 


8.877 


575 


7.324 


965 


10.133 



"Units: erg cm ^ s ^ A ^ at the stellar surface. 
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